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ABSTRACT 

Aims. We aim to investigate the effect of the escaping ionizing radiation on the color selection of high-redshift galaxies and identify 
candidate Lyman-continuum (LyC) emitters. 

Methods. We used the intergalactic medium (IGM) prescription of Inoue et al. (2014) and galaxy synthesis models of Bruzual & Char¬ 
iot (2003) to properly treat the ultraviolet stellar emission and the stochasticity of the intergalactic transmission and mean free path in 
the ionizing regime. Color tracks were computed by turning the escape fraction fesc of ionizing radiation on or off. 

Results. At variance with recent studies, a careful treatment of IGM transmission leads to no significant effects on the high-redshift 
broad-band color selection, even adopting the most extreme ionizing emission model (with an age of lMyr, zero dust, and metallicity 
Z/Z©=0.02). The decreasing mean free path of ionizing photons with increasing redshift further diminishes the contribution of the 
LyC to broad-band colors. We demonstrate that prominent LyC sources can be selected under suitable conditions by calculating 
the probability of a null escaping ionizing radiation. This was performed by running ad hoc Monte Carlo simulations anchored to 
the observed photometry, exploring the stochasticity of the IGM, and comparing the simulated and observed colors that encompass 
the Lyman edge. The method was applied to a sample of galaxies extracted from the GOODS-S field. A known LyC source at 
z = 3.795 was successfully recovered as a LyC-emitter candidate, and another convincing candidate at z = 3.212 is reported. A 
detailed analysis of the two sources (including their variability and morphology) suggests a possible mixture of stellar and non-stellar 
(AGN) contribution in the ultraviolet. 

Conclusions. The classical broad-band color selection of 2.5 < z < 4.5 galaxies does not prevent the inclusion of LyC emitters in the 
selected samples. High fesc in relatively bright galaxies (L > 0.1L*) could be favored by the presence of a faint AGN that is not easily 
detected at any wavelength. A hybrid stellar and non-stellar (AGN) ionizing emission could coexist in these systems and explain the 
tensions found among the UV excess and the stellar population synthesis models reported in literature. 

Key words. Galaxies: high - redshift; Galaxies: formation; Galaxies: distances and redshifts 


1. Introduction 


In the past two decades, photometric redshift and color selec¬ 
tion techniques have been widely used to identify high-redshift 
galaxies up to the reionization epoch, z > 6 ([Steidel et al. 1999 


Bunker et al. |2004HGiavalisco et al. |20Q4l|Vanzella et al. 


Castellano et al. |26lOHGrazian et al. |201 1 [|McLure et al. 2013 


2009 


Finkelste m et al. ||2012| |Bouwens et al. ||2014 ). The efficiency 

of the selection and the volume probed depend on the photo¬ 
metric system, the depth of the survey (Du ncan et al. ||2014|), 
the contribution of nebular emission (e.g., Schae rer & de Bar- 
ros |2010[ |Inoue |2011a| ) and the intrinsic star-formation history 
and duty cycle over whic h a galaxy becomes UV-bright and de¬ 
tectable in the ultraviolet ( [Wyithe et al. |2104| . Potential contam¬ 
inants include stellar sources, time-variable events like super- 
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novae, and spurious sources (e.g., Bouwens et al. 2014 


Inter¬ 
vening lower redshift objects that are not spatially resolved (or 
not de-blended) can also alter the observed magnitudes. There¬ 
fore, a well-defined cosmological volume (i.e., redshift inter¬ 
val) does not have a well-defined corresponding region in color 
space. 


The global shape of the spectral energy distribution and the 
discontinuities in the emerging electromagnetic spectrum have 
been used to estimate the redshift probability distribution func¬ 
tions (i.e., photometric redshifts). Typical observed breaks in 
the ultraviolet-optical regime are the Lyman edge at 912A, pro¬ 
duced by the photoelectric absorption of neutral hydrogen in 
the stars, the interstellar medium and the intergalactic medium, 
the Lya -break (1215.7A) caused by the redshift-dependent in¬ 
tergalactic transmission blueward of the Lya transition, and the 
Baimer break (A - 4000A), a proxy of the age of the stel¬ 
lar populations. The contribution of the nebular line emission 
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has also been demonstrated to be significant at high redshift 
by mimicking, attenuating, or boosting the expected underlying 
breaks (e.g., Baimer and/or Lycr-breaks, e.g .,|Schaerer & de Bar-| 
|ros ||2009| 2010; |de Barros et al.||2014[ |Vanzella et al. |2009[ ) or 
generating evident discontinuities (jumps) in the photometry, as 


measured in extreme emission-line galaxies (e.g., 

van der Wei 

et al. 20 lT} Amorin et al. ||2011; Maseda et al. 

2014) or in 


sources in which the signature of the Ha emission line or the 
[O hi] d4959,5007 + Hf3 emission structure have been detected 
in the Spitzer /IRAC channels at redshift z > 3 .8 and z > 6.5, 
respectively ( [Shim et al. |2011||Smit et al. |2013| ). Broad absorp¬ 
tion features like the 2175A bump are also expected to introduce 
dips in the photometry if a Milky Way-type extinction curve is 
present at high redshift (e.g.,|Capak et al |2011[|Kimm & Cen | 
|20l3) . 

Of these discontinuities, the Lyman edge (912A) in star¬ 
forming galaxies is still subject of investigation. Currently, a 
solid and direct measure of the stellar Lyman break at high red¬ 
shift (z > 2.5) has not been reported. This is relevant for char¬ 
acterizing sources responsible for hydrogen reionization and is 
still an active line of research at any redshift. A statistically 
significant sample of sources with a solid Lyman-continuum 
(LyC) detection and a comprehensive view on how it is spatially 
distributed relative to non-ionizing light are still lacking (e.g., 
Vanzella et al. ||2012[ |Nestor et al. ||2Q13[ |Iwata et al~| |2009 ; 


Inou e et al. 2Qllb||Siana et al. |2010| ). Qbservationally, this 


study is further complicated by the opacity of the intergalactic 
medium (IGM), which becomes significant with increasing red¬ 
shift. On the one hand, this opacity, coupled with the internal ab¬ 
sorption in the galaxy, allows us to select high-z galaxies through 
the Lyman-break technique mentioned above, on the other hand, 
it could prevent us from selecting sources with prominent es¬ 
caping ionizing flux. Recently, |Cooke et al. | ( |2014| ) concluded 
that a significant fraction of galaxies is missed in color-color se¬ 
lection diagrams if the ionizing radiation escapes. In practice, at 
redshift 2 < z < 5 the required drop in the Lyman-break selec¬ 
tion method could intrinsically introduce a bias and exclude LyC 
emitters from the selection. 


We have three main aims in the present work: (1) we address 
the above possible bias in the high-z selection by investigating 
the effect of the intergalactic transmission on the photometric se¬ 
lection of high-redshift star-forming galaxies when a leakage of 
ionizing radiation is present. To maximize the effect, we assume 
that all ionizing photons escape along the line of sight (fesc- 1). 
The nebular line emission and photometric scatter also perturb 
colors and alter the positions in the color-color plane and have 
been extensively explored and quantified in the past (e.g., pun- 
can et al.~||2Q14[ |Schaerer & de Barros ||2Q10l |Dahlen et al. 
2010]). These effects are also discussed. (2) We present a method 
for selecting LyC emitters by focusing on the detailed analysis 
of ultraviolet colors and deriving the probability of a null escape 
fraction of ionizing photons. This method allows us to identify 
potential LyC candidates without the use of a dedicated survey. 
It is presented and applied to a test sample of spectroscopically 
confirmed galaxies. (3) We test this selection method on a sam¬ 
ple of spectroscopically confirmed galaxies for which additional 
photometric and spectroscopic observations presented here al¬ 
low us to verify the results of our selection technique. 

In Sect. 2 the methodology is presented, and the results of 
the color-color selection techniques are discussed in Sect. 3. The 
method is applied to a sample of galaxies lying in the GOODS-S 
field, and the results are reported in Sect. 4. In Sect. 5 we discuss 
and summarize the results. 
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Errors are quoted at the ler confidence level unless otherwise 
stated. Magnitudes are given in the AB system (AB = 31.4 - 
2.5 log(/ v /nJy)). We assume a cosmology with D tot ,Q M ,Q A = 
1.0,0.3,0.7 and//o = 70 km s _1 Mpc -1 . 


2. Methodology 

We derived the synthetic colors of high-z galaxies by adopting a 
Monte Carlo prescription of the IGM transmission and including 
spectral templates to properly follow the ultraviolet stellar light 
emerging from galaxies (a schematic view is shown in Fig. [l}. 
In particular, we included two effects: (1) the treatment of the 
stochasticity of the IGM transmission both in the ionizing and 
non-ionizing regimes, in particular, the redshift evolution of the 
mean free path of ionizing photons is properly considered; (2) 
we used galaxy templates that cover the Lyman break and the 
ionizing emission, such that the shape of the Lyman continuum 
(A < 912A) is suitably represented and is self-consistent with 
the whole ultraviolet emission. We did not include dust attenua¬ 
tion for two reasons: first, although the dust extinction law would 
proceed monotonically in the ionizing regime and therefore we 
expect more attenuation in the LyC than at 1500A (e.g., Siana et 


[aT7| ( j2QQ7{ ), see their Figure 4), we prefer to avoid this extrapola¬ 
tion and assumed for now no dust effect on color-color selection. 
Second, when we exclude dust, we favor a stronger effect of ion¬ 
izing radiation on colors, therefore our approach is conservative. 


2.1. Intergalactic medium perscription 


|Inoue & Iwata [ |2008| (1108 hereafter) calculated the transmis¬ 
sion of the IGM by performing dedicated Monte Carlo (MC) 
simulations. These simulations have been used to correct for in¬ 
tergalactic attenuation and estimate the escape fraction of ion¬ 
izing radiation (e.g., Iwata et al.^] 2009[ | Vanzella et al.||2010b| 
2012). The MC simulation of 1108 have now been updated based 
on the new empirical function of intergalactic absorbers pre¬ 
sented in [Inoue et al. | ( |2014| ) (114, hereafter). We briefly recall 
the main steps. The simulations are based on an empirical dis¬ 
tribution function of intergalactic absorbers that simultaneously 
reproduces the observational statistics of the Lyman alpha for¬ 
est (LAF), Lyman limit systems (LLSs), and damped Lyman- 
alpha systems (DLAs). From this assumed distribution function, 
a large number of absorbers were generated along many lines of 
sight. The probability of encountering an absorber was assumed 
to follow a Poisson distribution, and for each one the column 
density and Doppler parameter were extracted randomly from 
their (empirical) probability distribution functions. 

The simulations are therefore anchored to the observed dis¬ 
tribution of th e number density of abs orbers at vari ous column 
densities (e.g., Prochaska et al. 2009] 2010, 2014; |Songaila~&| 
|Cowie ^OlOUQ’Meara et al. |2013|), and verified to be consistent 
with the mean Lyg transmission ( [Becker et al. ||2013[ [Faucher- 

Gig uere et al. ||2008|) and with the mean free path of ionizing 
photons (Worseck et al. 201 A) in the redshift interval 1-6 (see 
114 for more details). We also verified that the mean and scat¬ 
ter at A > 912A was in line with that observed in QSOs stud¬ 
ies (e.g., Faucher-Giguere et al. |||2008] [Rollinde et al. I|20l3| ). 
The distributions of the median transmissions with the central 
68% intervals computed in the rest-frame ranges 880-910A and 
1070 - 1170A are shown in Fig. [^J in which the asymmetric dis¬ 
tribution is evident in both cases with a tail toward low transmis¬ 
sions due to the presence of DLA and LL systems (see 114 and 
1108 Fig. 7 for details). As discussed in 1108, the absorption of the 
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LyC is mainly caused by LLSs and DLAs with Nhi > 10 11 cm~ 2 . 
This implies that LyC absorption is very stochastic because it is 
related to the probability of intercepting dense absorbers (e.g., 
LLSs). The simulations (which are based on the new absorber 
function by 114) used in this work consist of 10000 random lines 
of sight at a given redshift, in the range 0-7 with dz = 0.1, 
each transmission spans the wavelength (rest-frame) interval 
600-1300A with dA= 1A (an example of a single l.o.s. is shown 
Fig. [T]). Redward of the Lya transition (1216A), the IGM trans¬ 
mission is assumed to be 1.0. These simulations are therefore 
ideal to account properly for the effect of the IGM stochasticity 
on color selection. 


2.2. Galaxy template spectra 

The ultraviolet rest-frame spectra of star-forming galaxies are 
well reproduced with spectral synthesis models (e.g., |Ilbert et al.| 
2009}[Heckman et al. [201 1 [Leitherer et al. |2014 }. We a dopted 


galaxy template spectra of |Bruzual & Chariot ( |2003| )(BC03, 
hereafter), which are widely used in the literature and clearly de¬ 
scribe the typical ultraviolet and ionizing emission arising from 
star-forming galaxies. At the level of accuracy required here, 
more sophisticated models of stellar evolution do not introduce 
significant differences in the results (e.g., |Levesque et al. |2012| ). 
We must rely on theoretical models because of the lack of ob¬ 
served spectra of confirmed galaxy with Lyman-continuum leak¬ 
age below the Lyman limit. 

The shape and intensity of the ionizing radiation is linked 
to the nature of the stellar populations involved, such that 
for standard stellar populations the luminosity density ratio 
L(1500A)/L(T) increases as wavelength decreases (A < 912A, 
e.g., [Inoue et al. ||2005| |Siana et aL||2007| ). An extreme case is 
shown in Fig. [jj in which a galaxy template placed at z = 3.0 
with zero dust extinction, E(B-V)=0, low metallicity Z/Z© = 
0.02 and young, lMyr old, is superimposed to the photometric 
filters used in this work. These parameters correspond to an ul¬ 
traviolet slope of (3jjv - -3.3, which is quite extreme and very 
rarely observed (e.g., [Vanzella etal. |2014a| . The ultraviolet den¬ 
sity ratio for this template is F y (1500)/F y (900) = 1.7, extremely 
blue if compared with typical star-forming galaxies at the same 
redshift (F y (1500)/F y (900) ^ 6 , |Siana et al.~| ( 20 T 0 )). However, 
we used this extreme case to maximize the effect of the ionizing 
emission on colors when diX\fesc=\ is assumed. Even more ex¬ 
treme ionizing emission can be realized as in the Lyman-bump 
scenario described by Inoue ( |2010| ). Here we limit our analysis 
to the standard BC models. In the same figure a random inter- 
galactic transmission (one out of 10000 at z = 3) is also shown 
(magenta solid line). In this case a relatively high IGM trans¬ 
parency is present down to - 600A rest-frame. A galaxy tem¬ 
plate with a redder and more typical UV slope (fiuv - -2.1) has 
also been adopted (E(B-V)=0, Z/Z© = 0.2, age=lGyr). 


3. Results of the color selection 

U n — G,G — *R ( [Steidel et al. |2003[ ), $435 — V 606 ? and V 6 O 6 — his 
(HST/ACS) colors were calculated at each redshift (2.7 - 4.8, 
dz = 0.1) by convolving 10000 IGM l.o.s. with filters and BC 
galaxy templates. The escape fraction of ionizing radiation has 
been turned off or on by imposing zero signal blueward of the 
Lyman edge (fesc=0 ) or by retaining the shape of the original 
BC template (fesc=l). The emerging ionizing flux is consistent 
with the rest of the ultraviolet emission, meaning that it is not 
artificially added. 



Fig. 1 . Examples of the galaxy template (F v , black line) placed at red¬ 
shift z = 3.0, intergalactic transmission (one l.o.s. out of 10000) and 
the filters used in this work. The Lyman continuum is indicated (shaded 
region), and in this example it affects half of the U n - band. The BC tem¬ 
plate shown is the most extreme used in this work and maximizes the 
ionizing emission (it corresponds to an ultraviolet slope of f3 = -3.3). 



redshift redshift 

Fig. 2. IGM transmissions as a function of redshift calculated in the 
wavelength intervals 880-910A (left) and 1070-1170A (right). At each 
redshift, 10000 line of sights have been computed (black points), and 
the medians with the central 68% intervals are also shown (red filled 
circles with errorbars). 


Figure [3] shows the results in the U n GR (U n - G vs. G - R) 
and BVI (#435 - V 606 vs. V^o6 ~ his) color-color planes typi¬ 
cally used to select galaxies in the redshift range 2.5 < z < 3.5 
(U n - band dropouts, Steidel et al. ||2003|) and 3.4 < z < 4.5 
(# 435 -band dropouts,"]Giavalisco et al. |2004| ). To set limits in 
the color diagrams, a one-sigma limiting magnitude of 30 (in all 
bands) was assumed and the template spectrum was normalized 
to £850 = 25.0. The saturation in the upper part of the color-color 
diagrams visible as a decreasing linear envelope is due to this as¬ 
sumption. In this exercise we did not include photometric scatter 
(but see the discussion below). While the dust absorbs the UV 
radiation, both metallicity and age anticorrelate with the produc¬ 
tion of ionizing radiation, such that the younger the stellar popu¬ 
lation and/or the lower the metallicity, the higher the production 
of ionizing photons (e.g., |Nakajima & Ouchi |20f4| . 

Two BC templates were adopted and are shown in Fig.[3j the 
extreme case already mentioned above (J3 = -3.3) aimed to max¬ 
imize the effect of LyC on colors, and another more typical case 
(even though it is still relatively blue, = -2.1). We indicate in 
the same figure with an arrow (direction and magnitude) the ef¬ 
fect that E(B - V) = 0.2 would have on colors by monotonically 
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Fig. 3. Top panels: The U n GR color planes with the typical selection regions for galaxies (solid blue line) in the redshift range 2.5 < z < 3.5 are 
shown for moderately blue (left) and extremely blue (right) galaxy templates. The effect of dust extinction and Lya emission is shown with red 
arrows. The green clouds and black stripes correspond tofesc=0 andfesc=l , respectively, and are computed with redshift steps of 0.1 indicated 
at the left side. Magenta and cyan solid lines connect the barycenter of each cloud/stripe and represent the average color tracks (see text). In the 
top-left panel the (solid) black color track with filled (red) squares denotes the case in which the IGM is absent from the ionizing part of the 
spectrum, i.e., the transmission is equal to 1. The rightmost filled square symbol represents z = 4.5; squares have been plotted with dz = 0.1. 
The color track resembles the behavior of the flat models of Cooke et al. (2014) (see Sect. 3.1). Magenta dashed lines underline the BM and BX 
selection regions. Bottom panels: The same as top panels, but in this case, the BVI (5 435 , V 60 6, ins) color selections of galaxies expected to lie 
in the redshift range 3.4 < z < 4.5 are shown (solid blue lines). In the bottom-left panel the star indicates the position of the Lyman-continuum 
emitters Ionl at z - 3.795 (see text). In the bottom right panel the vertical orange dotted line highlights the bending of the black stripe at z = 4.5 
when ionizing emission also affects the V 60 6-band, visible as a slightly bluer V 60 6 - 1775 color for the higher IGM transmission cases (see text). 
The effect of the LyC leakage on colors decreases with increasing redshift (passing from U -band to 5-band dropout selection methods) as a result 
of redshift evolution of the mean free path of ionizing photons. In all quadrants the triplets of red arrows mark the direction and magnitude of 
the color modification as a result of the Lya emission, nebular continuum (affecting only thefesc=0 case, green points), and the dust attenuation, 
E(B-V). 


extrapolating the law of Calzet tTet al. | ( |2000| ) (see also, Siana 


|et al.||2007] ). The green and black points shown in Fig. |3] were 
calculated adopting fesc=0 and 1 , and are named green clouds 
or black stripes in the following. 

[fesc=0]: Each green cloud shown in Fig. [3] is made of 
10000 points (that is, 10000 l.o.s.), and their behavior repre¬ 
sents the color regions with fesc=0 at various redshifts. The 
redshift intervals explored in the U n GR and BVI diagrams are 


2.7 < z < 4.1 and 3.4 < z < 4.8 (with dz = 0.1), respectively. 
With increasing redshift, the Lyman edge and the null flux at 
A < 912A enters the U n - band, producing a rapid drop (# 435 - 
band in the BVI case). The spread in the horizontal and vertical 
directions reflects the IGM scatter in the non -ionizing domain, 
A > 912A. In the horizontal axis it only affects the G-band (V- 
band in the BVI case), while in the vertical axis both U n and G 
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bands are affected. This scatter is in line with that observed in 
QSOs studies (e.g., Faucher-Giguere et al. 2008t [Rollinde et al 


2013; 114). Focusing on the U n GR plane, the clear ( U n -G ) color 


cut in the bottom envelope of the green clouds is therefore due 
to a sharp break at 912A (fesc=0 ), and the scatter in the (G-R) 
reflects the scatter of the IGM at (A > 912A). The barycenters 
of the green clouds (connected with the magenta line in Fig. [3]) 
resemble the typical color tracks reported in the literature, ac¬ 
cording to the adopted galaxy template and dust attenuation. 


\fesc=l]: More relevant in this work is the case withfesc=l. 
The stochasticity of the IGM transmission in the LyC region and 
its effects on colors is appreciable. For example, the tail of trans¬ 
missions toward high values (transparent l.o.s.) is clearly visible 
in the U n GR diagram as a percolation of points toward relatively 
blue ( U n -G ) values (vertical black stripes). This is more evident 
in the extreme case (Puv = -3.3), in which the ionizing emission 
is maximized in the set of models we used. We note that the me¬ 
dian color tracks for fesc-O/l can differ in the U n -G color in the 
U n GR plane by a negligible quantity up to one magnitude in the 
P = -2.1 case. It is even higher for the most extreme template 
((3 = -3.3) because the ionizing flux is higher. This difference 
decreases as redshift increases (because of the decreasing mean 
free path of ionizing photons with redshift), as shown in the BVi 
z - 4 plane (see Fig. [3]). However, and most importantly, the es¬ 
caping ionizing radiation does not significantly move the points 
outside the selection regions. In particular in the ft = -2.1 case, 
there are no solutions outside the selection region as a result of 
the ionizing emission. The solutions escaping the region in the 
low-redshift tail, z < 2.8 (z < 3.5 in the BVI case) are due to 
the intrinsic blue slope, regardless of fesc. In the more extreme 
case,/? = -3.3, only a marginal fraction leaves the selection box 
( U n GR ) at z > 3.7, a redshift range that is also covered by the 
other higher-z BVI color selection, however. In the U n GR and 
BVi diagrams the fraction is always < 10% in the p = -3.3 BC 
template. However, the selection boxes were designed to capture 
star-forming galaxies with typical properties, therefore it is not 
surprising that the p = -3.3 case does not exactly match the 
redshift range, especially the lower redshift edge of the selection 
volume. We recall that the two templates adopted here are dust- 
free and tend to maximize the contribution of ionizing emission 
to the colors. 


The effect of an escaping ionizing radiation (fesc= 1) de¬ 
creases as redshift increases, mirroring the redshift evolution of 
the mean free path of ionizing photons. For this reason, the same 
black stripes in the BVI plane (redshift 3.4 - 4.5) are less promi¬ 
nent than in the lower z U n GR plane. If the m.f.p. and its red¬ 
shift evolution are neglected from the calculation, the dropout 
signature is strongly attenuated and color tracks escape from the 
selection region (see next section). 


A second-order effect of fesc- lisa small bending of the 
stripes in the case of high IGM transmissions if the Lyman edge 
affects both colors in the given color-color diagram ( U n GR or 
BVi). An example is indicated with a vertical orange dashed line 
(Rg. a bottom right panel) superimposed to the black stripe 
at z = 4.5. A bending in the bottom envelope of the stripe at 
(#435 ~ V 506 ) - 2.0 and (V 606 - his) - 0.7 is appreciable. This 
is because the LyC also boosts the V 606 -band, producing slightly 
bluer values (V606 - z'775). However, this is a second-order effect 
that can be easily erased by photometric noise. 


3.1. Marginal effect of the escaping LyC on the broad-band 
color selection 


As shown above, the effect of an escaping ionizing radiation on 
the color selection is minimal. Conversely, [Cooke et al. ( |2014| 
found a significant contribution and calculated color tracks that 
lie outside the selection regions when the escaping ionizing ra¬ 
diation is present. They performed spectrophotometry on com¬ 
posite spectra based on z ~ 3 LBGs (e.g., |Shapley et al. |2003| ) 
and applied various levels of escaping ionizing flux down to a 
Amin, and assuming various /1500 / f Ly c ratios. Two models have 
been adopted: 1) the case in which A M in is not specified (called 
flat-model) and 2) the case in which A M in = 800A rest-frame 
(called step-model). While the adopted non-ionizing spectrum 
is (by construction) close to the average of the observed pop¬ 
ulation, the assumptions on the Lyman continuum are too sim¬ 
plistic and lead to unrealistic color tracks, especially in the flat- 
model. The flat-model neglects the m.f.p. and its redshift evo¬ 
lution, therefore the dropout signature is significantly reduced. 
The color tracks plotted over the U n GR plane in Cooke et al. 
(2014) (see, e.g., their Fig. 3 and 4) show colors outside the se¬ 
lection region, reaching the extreme case at z = 4.5 in which the 
U n - G is close to zero, that is, with increasing redshift, the LyC 
enters the G-band and still survives in the U n band. At z = 4.5 
the G„-band probes A < 730A rest-frame, while the m.f.p. at 
this redshift is 22Mpc (Worseck et al. |2014| ), and corresponds to 
880A rest-frame, that is, on average, an ionizing photon survives 
on a much shorter path than is allowed in the flat-model. The 
top-left panel of Fig. [3] shows the color track of the BC template 
with p = -2.1, fesc=l and assuming 100% IGM transmission 
blueward of the Lyman edge (A < 912A, solid black line and red 
squares). The behavior of the flat-model is well reproduced (up¬ 
turned shape) with the dropout in the G„-band that never occurs 
up to z = 4.5 (as in Cooke et al.). As discussed above, neglect¬ 
ing the m.f.p. of ionizing radiation in the intergalactic medium 
leads to a strong attenuation of the observed drop. The compar¬ 
ison between the cyan and red tracks in the same figure shows 
the strong effect that is achieved when the IGM is neglected in 
the ionizing domain. The step-model is more realistic than the 
flat-model and produces color tracks typical of what has been 
reported in literature (e.g., Steidel et al.), although, again, it does 
not include the redshift evolution of the m.f.p. of ionizing pho¬ 
tons. The m.f.p. at z = 3.0/3.5/4.0 is - 100/50/33Mpc proper 
( jFumagalli et al. [2013[|Worseck et al. |2014[ ), which corresponds 


to Amin ~ 820/858/870A rest-frame. Therefore the assumption 
of Amin = 800A is valid only for a narrow redshift range around 
z < 3 and tends to underestimate the growth-rate of the U n - G 
color with increasing redshift (z > 3). 

We here properly treat the aforementioned limitations, that 
is, the shape of the LyC flux (modeled with galaxy BC- 
templates) and the m.f.p. of ionizing photons (statistically de¬ 
scribed by the 108,114 prescription) are both included in the MC 
simulations. Our results show that there is no significant effect 
of ionizing leakage on the broad-band color selection. However, 
as discussed in Sect. 4, a leakage of ionizing radiation is in prin¬ 
ciple identifiable under suitable conditions. 

Another effect that can decrease the photometric breaks or 
change the UV slope in galaxies and move points toward the 
outer part of the selection regions is the (photometric) contami¬ 
nation by lower redshift superimposed sources, which is not neg¬ 
ligible in the case o f deep and/or low spatial resolution surveys 
(e.g., see Vanzella et al. | ( |2Q10a| ) for a statistical analysis). Taken 
as a whole, it is worth noting that this type of contamination 
is different from the effect produced by a pure escaping ioniz- 
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ing radiation. The latter is modulated by the redshift evolution 
of the IGM opacity (m.f.p.) and therefore limited in its global 
contribution to the band. Conversely, the contribution of a lower 
redshift source mimicking the LyC emission typically affects the 
entire band(s) blueward of the Lyman edge, which significantly 
reduces the observed drop or causes the background higher-z 
galaxy to not drop at all, similarly to what was calculated in the 
flat-model of Cooke et al. 


3.2. Photometric noise, nebular continuum, and emission 
lines 


The effect of photometric noise on color selections has been ex¬ 
tensively explored in the literature (e.g., |Dahlen et al.~||2010 , 
|Duncan et al. |[20l4| ). The photometric scatter depends on the 
signal-to-noise ratio (S/N) in the bands, which in turn is related 
to the depth of the survey and/or to the magnitude of the sources. 
In the present case - at fixed conditions - the fesc=0 produces 
fainter U n -band (or 5435 -band) magnitudes than thefesc=l (by 
definition), therefore we expect that thefesc=l case is slightly 
less scattered than the fesc=0. It is beyond the scope of the 
present work to address this_question or replicate the analysis 
performed, for instance, by [Dun can et al. ( j2014| ). We only recall 
that depending on the adopted depth of the survey and magnitude 
of the sources, the photometric scatter can move poi nts o utside 
the selection boxes, and vice versa. Dun can et al. | ( 12014] ) found 
that most of the cases that escape the selection regions are photo¬ 
metric limits, meaning that they are compatible with the high-z 
locus in the color space (with virtually infinite photometric ac¬ 
curacy, they would be included in the high-z regions). Clearly, 
the high-z selection is futile above a certain level of noise. 


The Lya line in emission can also affect the observed colors. 
Since the contribution to the broad-band is proportional to the 
observed equivalent width of the line, this effect is more signifi¬ 
cant with increasing redshift. Given the redshift interval explored 
here (2.5 < z < 4.5), the contribution of the Lya emission is rela¬ 
tively small and anyway favors the color-color selection criteria, 
that is, it tends to move points toward the selection region. For 
th eU n GR color plane the Lya line always lies in the G-band and 
therefore any emission moves points to the upper left part of the 
diagram. As an example, the equivalent width of 65 A rest-frame 
(corresponding to the highest theoretically expected value in the 
P = -2.1 case adopted here and assuming case B recombination) 
introduces a differential magnitude dm - 0.15 in the ( U n - G) 
and ( G-R ) colors, in the redder and bluer directions, respectively 
(the effect is shown as a vector in Fig. [3]). The same occurs for 
BVI . Clearly, if the Lya line enters the U n - band (or #435 -band 
in the BVI case), the effect is the opposite: the line decreases 
the observed break. However, this is not the case here. In the ex¬ 
treme case adopted here (with ft ^ -3.3), the effect of a Lya line 
with rest-frame equivalent width of 375A is more significant, as 
expected. The nebular continuum can also redden the UV slope 
(e.g.,[Castel lano et al. |2014 ). It is maximized in the/Lsz^O case, 
in which the ionizing radiation concurs to boost flux in the non¬ 
ionizing UV range (Schaerer 2002; |Inoue et al. |2011b| ). In the 
moderate case (J3 = -2.1), the nebular continuum effect is neg¬ 
ligible, producing a A/3 = +0.07, while in the ft = -3.3 case a 
A p = +0.60 is obtained, clearly in line with the extreme physi¬ 
cal conditions (the effect on colors is reported with an arrow in 
Kg# 
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4. Photometric selection of LyC candidates 

Ionizing sources have been searched for by directly examining 
the Lyman continuum emission at 2 < z < 4.5 with deep spec¬ 
troscopy ( [Steidel et al. |2001[|Giallongo et al. |2002| |Shapley et 


al. |2006 ), narrow-band imaging (e.g., Iwata et al. 2009, Nestor 

et al. 12013} |Mostardi et al. ||2013|), and deep inte rmediate and 


broad-band imaging (e.g., Femandez-Soto et al. |2003[|Siana et 

|al.|2007| 2010; |Vanzella et al.|2010bl 2012; |Boutsia et al. |2011i 
Grazian et al. 2015). Deep slitless spectroscopy has also pro¬ 
vided deep limits at z ~ 1 ( Bridge et al. |2010| . 

As discussed in the previous section, the leakage of stellar 
ionizing radiation (fesc= 1) convolved with the IGM does not 
compromise the color selection of high-z galaxies. Despite this, 
a statistically significant sample of LyC emitters (fesc>0 ) at high 
redshift is still lacking (e.g., see |Vanzella et al. |2012j ). One rea¬ 
son could be that a significant number of ionizing sources exist 
only in a much fainter and still unexplored luminosity domain 
(e.g., [Wise et al. ||2014[|Kimm & Cen ||2014 ), while relatively 
bright (L > 0.1L*) ionizing sources are rare. The contribution to 
the ionizing background therefore could be negligible if we fo¬ 
cus on the bright part of the luminosity function, but identifying 
them is still interesting since the physical mechanisms that al¬ 
low ionizing photons to escape are not fully understood and are 
probably linked to efficient feedback processes (e.g., AGN). 

In the following we describe a method for identifying LyC 
candidates and apply this to a sample of spectroscopically con¬ 
firmed galaxies in the GOODS-S field. We defer a systematic 
study on a larger spectroscopic redshift sample to a forthcoming 
work. 


4.1. Searching for anomalous UV colors as a signature of an 
escaping ionizing radiation 

As shown above, an optimal characterization of the stochasticity 
of the IGM and a suitable modeling of the ultraviolet emission 
of star-forming galaxies allows us to investigate the photomet¬ 
ric contribution of an escaping LyC radiation in more detail. In 
particular, if the photometric accuracy is high enough, it is possi¬ 
ble to explore the position of the source in the color-color plane, 
focusing on the color(s) encompassing the Lyman edge, and de¬ 
termine whether it significantly departs from thefesc=0 assump¬ 
tion. 

To this aim, we built a set of templates for which we adopted 
different values of the ages, metallicities, and dust attenuation. 
Here the dust attenuation is inserted because we compare sim¬ 
ulated with observed colors (discussed in the next sections). 
The following grid was adopted: Z/Z Q = 1,0.2, and 0.02; E(B- 
V)=0.0,0.03,0.06,0.1,0.2,0.3,0.4,0.5, and 0.6; age=lMyr up to 
2Gyr with steps of 0.05 in logarithmic scale. Both constant 
and exponentially declining (with r = 0.6Gyr) star formation 
histories were assumed. The Salpeter initial mass function and 
Calzetti attenuation law extrapolated down to the LyC were 
adopted (e.g., |Siana et al.|2007| ). The final sample uniformly cov¬ 
ers the age, metallicity, and dust attenuation parameter space. 

The basic idea is the following. Given a galaxy with known 
spectroscopic redshift (zspecj] the UV spectral slope (p o bs) 
with its one-sigma uncertainty (erg) is measured by computing 
a multi-band fitting as described in Castell ano et al. | ( |2012| ). 
The uncertainty on the estimated slope, crp, contains the pho¬ 
tometric uncertainties. The same procedure is performed on BC- 


1 We have shown that the photometric selection does not preclude the 
inclusion of LyC emitters. 
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templates to derive the UV slope ( p mo dei ) at the same redshift and 
using the same bands. Subsequently, a subsample of galaxy tem¬ 
plates with f3 mo dei compatible with the observed one is selected 
such that p obs - crp < p mo dei < Pobs + erg. To be conservative, we 
furthermore select in this sample the template that maximizes the 
ionizing output, by measuring the Ll500/L LyC flux density ratio. 
This also means that we select the template that maximizes the 
contribution of the non-ionizing ultraviolet emission redward of 
the Lyman edge (e.g., > 912A). This part still affects the band 
even iffesc=0. In this way, we build a galaxy template anchored 
to the observed UV slope and with the highest possible (mod¬ 
eled) ultraviolet emission. Then we convolve this template with 
10000 IGM transmissions at the closest available redshift in the 
IGM redshift grid, zigm (always zspec-ziGM < 0.05) and derive 
magnitudes and colors, turning the escape fraction on and off, 
and normalizing the template to the observed zsso magnitude. In 
particular, we focus on the realization that assume fesc=0 and 
compare them with the observed photometry, with the aim to 
identify a possible excess of the flux in the ionizing domain. 

The hypothesis is that we have no direct coverage just blue- 
ward of the Lyman edge (912A), therefore we have no direct 
probe of the LyC region, as is the case in general for non - 
dedicated surveys. We wish to infer a signature of the ionizing 
radiation by focusing on the band in which the LyC is partly 
covered. While dedicated surveys based on narrow-band imag¬ 
ing and deep spectroscopy are more efficient in capturing the 
possible escaping ionizing flux just blueward of the Lyman limit 
(although they are time-consuming and/or limited to narrow red¬ 
shift intervals), the method described here complements them 
and allows us to explore possible LyC leakage in a continuous 
redshift range, exploiting existing deep multifrequency surveys. 
The mean free path of ionizing photons is ~ lOOMpc physical at 
redshift 3 ( [Fumagalli et al. |20l3[ [Worseck et al. |2014| ), which 
corresponds to ~ 80A rest-frame blueward of the Lyman edge 
at the same redshift, that is, ~ 320A in the observed-frame. If 
compared to the typical width of broad-band filters, this is not 
negligible (but not dramatic for the high-z selection functions, 
Sect. 3). 

For example, the (U)B ^-band starts to include the LyC do¬ 
main at z >(2.30)3.05, and at z =(3.0)3. 8 , half of the ( U)B^s- 
band probes ionizing radiation. The distributions of the simu¬ 
lated colors and magnitudes derived assuming fesc- 0 and con¬ 
taining the Lyman edge (912A) are then compared with the 
observations. In the case in which the # 435 -band contains the 
Lyman limit, the probability that at the given observed color 
(B - V) obs , UV slope p and spectroscopic redshift (with their 
uncertainties) is associated with anfesc=0 ( P(fesc=0)) is calcu¬ 
lated as follows: 


P(fesc— 0)\[z spec , P obs ±lo- p , (B-V) obs +lo- (B -v)] - N[color]/N[totaI], (1) 

where N[ tota i] = 10000 and N[ co i or ] is the number of real¬ 
izations with/<?sc =0 that reproduces the observed color, that is, 
C B - V) simu i <{B- V) obs + 1 cr (B - V ), calculated adopting the BC 
galaxy template with the strongest LyC emission of those that 
are compatible with the observed ultraviolet slope p obs (±1 erg). 
Therefore, given the observed quantities (color, redshift, and UV 
slope), a proper treatment of the IGM transmission (114) and UV 
emission (BC-models), the probability of a null fesc can be es¬ 
timated from MC simulations. In particular, the LyC candidates 
are those with the lower conditional probability P(fesc=0). 

We applied this procedure to a sample of galaxies in the 
GOODS-S field with known spectroscopic redshift ( [Vanzella °et 


|al. ||2008[ 2009; [Balestra et al~||2010| ). In this first test-run we 
focused on the 2.9 < z < 4 redshift interval. The lower redshift 
limit was chosen such that more than half of the (ground-based) 
U-band probes the LyC. ^ The upper limit is mainly given by 
the short mean free path of ionizing photons at high redshift, 
z > 4 (m.f.p. < 35Mpc proper), which makes the contribution 
to broad-band filters negligible on average. For each source the 
P(fesc=0) was calculated using the (U — B) and (B-V) colors as 
described above. At z > 3.4 the Lyman limit moves redward of 
the red edge of the U -band filter, which means that the U -band 
is a direct probe of the LyC. These cases have previously been 
explored in detail in Vanzella et al. (2010b), where only one con¬ 
vincing LyC emitter (out of 102 LBGs) was identified and named 
lonl (see below). Several other U-band detections at z > 3.4 
have been identified in the same GOODS-S field, but a care¬ 
ful investigation that exploited the optical and near-infrared high 
spatial resolution imaging (GOODS+CANDELS) has revealed 
a probable contamination by foreground sources that mimics the 
LyC (each of them analyzed in Vanzella et al. |2012] ). The frac¬ 
tion of contaminants is also compatible with that expected from 


statistical arguments (see Vanzefia et al. 2010a| 


In this application we furthermore selected sources with a 
S/N > 5 in all the ultraviolet bands and with secure spec¬ 
troscopic redshift measurements. We obtained a sample of 32 
galaxies in the range 2.9 < z < 4, magnitude 24 < zsso < 25.5, 
and ultraviolet colors redward of the Lya typical of LBG popu¬ 
lations (e.g., Vanz ella et al. |2009] ). The probabilities we derived 
(. P(fesc=0)) range between 0 and 100 %, with some cases of fore¬ 
ground contamination (discussed in Vanzella et_al. | |2012[ ) and 
known X-ray detected AGNs. In the following we focus on the 
sources with the lowest probability P(fesc=0)< 5% that are iso¬ 
lated, that is, those without evident contamination from nearby 
lower z objects, deferring a systematic application on a larger 
spectroscopic sample and a detailed comparison of LyC emitter 
candidates and non-candidates to a future work. Here we aim to 
demonstrate the feasibility of the method. In this way, two out of 
32 sources show P(fesc=0)< 5%: one is the already known lonl 
source at z = 3.795 (see next section and Vanzell a~et al. |2Q12| ), 
the other is a new intriguing candidate ionizer (see Sect. 4.3). 
These two objects are investigated in more detail in the follow¬ 
ing. 


4.2. LyC emitter: lonl 


The procedure described above recovers lonl as a LyC emitter 
candidate (CANDELS ID 23836, |Guo et al. |2013| ), even exclud¬ 
ing the direct information we have about the LyC detection (i.e., 
the U-band detectior0at A < 830A rest-frame, Fig. [3]). 

We recall that, currently, the global spectral energy distri¬ 
bution of this source is compatible with a star-forming galaxy: 
there is no X-ray detection in the 4Ms Chandra observations, 
and no emission lines (e.g., Lya, NV, CIV) have been measured 
in the spectrum ( Vanzella et al! |2010b[ 2012). It is a compact 
source, but is marginally resolved in the ACS images (A ~ 1500 
rest-frame) with a R e ff - 250 pc (see Sect. 4.6). 

At the redshift of lonl more than half of the #435 band probes 
the ionizing part of the spectrum (< 912A) (see top-left panel 


2 Bluer bands like the ultraviolet channels provided by the 
HST/WFC3, e.g., F225W, F275W, and F336W, will allow us to investi¬ 
gate lower redshift regimes. 

3 Here the U -band refers to the VIMOS U -band, Nonino et al 


(2009! 


whic h is an intermediate-band filter, FWHM=350A (see also Vanzella 
et al.|2010b| ). 
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Fig. 4. Summary of the lonl source. In the top-left panel the HST/ACS, 
VLT /U, and 5-band filters are shown, superimposed to a BC galaxy 
template placed at the same redshift of lonl. In particular, the solid blue 
and red dotted lines represent the HST/5 435 and VLT/5 band filters. The 
shaded transparent blue region marks the LyC region. More than half of 
the 5 435 and VIMOS/B bands probe the ionizing domain, while the U - 
band probes A < 83OA rest-frame. In the top right panel the simulated 
distribution of the (5 435 - V 6 oe) colors are calculated assuming fesc-0/l 
and compared to the observed color of the source lonl . The probability 
of a null fesc is very low. Similarly, the simulated U -band magnitude 
reproduces the observed one only marginally (middle right panel). The 
bottom panel shows the thumbnails for lonl . The blue circle in the 
#435 -band has a diameter of 1". 


of Fig. [4]). The UV slope measured with the five-band fitting 
procedure gives p o b s = -2.01 ± 0.10. The observed color is 
(#435 ~ V( 506 ) = 1-22 with an error cr^-v) = 0.18, the zgso-band 
magnitude is zsso = 24.83 + 0.03. We derived the probability 
P(fesc=0) by calculating the expected (#435 - V606 )simui color 
including the 114 IGM prescription and adopting the galaxy 
(BC) model anchored to the observed UV slope, as discussed 
in the previous section. The distribution of the expected color 
(#435 - V606) simul for the lonl source is shown in Fig. [4] and 
is compared to the observed color. Interestingly, although this 
source was selected as a #435 -band dropout (indicated with a star 
in Fig. [3] bottom-left panel), the observed (#435 - V606) color is 
reproduced only marginally iffesc=0 , that is, P(fesc=0) - 4.7% 
(shaded region of the fesc- 0 histogram shown in Fig. 0 up¬ 
per right panel). With such a low probability it is selected as a 
LyC emitter candidate. Conversely, the color is better reproduced 
with the distribution that assumes fesc- 1 (see the red curve in 
the same figure). We recall that the adopted template is the most 
ionizing while reproducing the observed UV slope. If we run the 
same procedure on all the templates that satisfy the condition 
Pobs - 0-/3 < Pmodei < Pobs + ct/3 (239 BC templates) regardless 
of any constraint on the ionizing emission, 2.39 x 10 6 color esti¬ 
mates (= N[total]) are derived. In this case, P(fesc=0) is 0.4%. 

A similar analysis can be performed on the observed mag¬ 
nitudes that probe the LyC domain. As described above, the 
galaxy template was normalized to the observed zs 50 • While the 
expected magnitude in the #435 -band better matches thefesc=l 
distribution (although its error still prevents us from clearly dis¬ 
tinguishing between fesc- 1/0), the observed U-band magnitude 
of 27.92 ± 0.16 - that probes the A < 830A rest-frame - is 
more difficult to reproduce fesc- 1 and given the observed UV 
slope, even in the more transparent IGM l.o.s. andfesc=l , see 
Fig. [4](the case withfesc=0 produces zero flux in the U-band by 
definition). The measured U-band flux can be better explained 


Article number, page 8 of 13 




27.0 27.5 28.0 

UviMOS 


•i.vS 1 

VLT/Uvimos B435 Veoe 


Fig. 5. Summary of the Ion2 source, see description of Fig. [5] The slit 
orientation is shown in the / 775 image with the wavelength dispersion 
along the vertical direction. 


with bluer and more ionizing BC templates (e.g., p < -2.5), 
which contradict the observed p = -2.01 ± 0.10, however. This 
could suggest a non-stellar radiation that progressively domi¬ 
nates at shorter wavelengths (especially in the ionizing domain). 
In particular, if we adopt a simple underlying power-law tem¬ 
plate without any break at the Lyman edge (AGN-like), F A = Af 
with p = - 2 , we obtain even bluer (#435 - V 606 ) colors and 
brighter U -band magnitudes than observed. Obviously, pump¬ 
ing flux into the ionizing regime allows us to recover the U- 
band magnitude as well. We recall that this object is currently 
not detected in X-ray (the 4Ms Chandra observations provide a 
lcr limit of Lx < 3.5 x 10 42 erg s~ l ) or in the 24 pm Spitzer/MIPS 
observations ( [Vanzella et al. |2012] ), and it is spatially resolved in 
the HST/ACS images (see Sect. 4.6). A possibility is that a mix 
of stellar and (faint) AGN emission contributes to the ultraviolet 
spectrum. More in general, it could be that relatively bright ion¬ 
izing sources host a faint AGN that favors the escaping ionizing 
radiation by decreasing the column density of neutral hydrogen 
gas along the fine of sight, causing fesc to be higher than zero, as 
in the lonl case (i.e., Nhi < 10 18 cm -2 ). This would imply that 
faint AGNs can still provide the necessary feedback to increase 
the transparency toward certain view angles (e.g., Giaflongoetj 
al. |2012| ). The ionizing photons produced by young stars in the 


compact central region (revealed in the UV at A > 912A as spa¬ 
tially resolved emission, see Sect. 4.6) could share a similar (or 
the same) physical path of non-stellar ionizing emission (AGN) 
and escape from the interstellar and the circum-galactic media. 
In this scenario, the final escaping ionizing radiation would have 
a stellar and non-stellar origin, with a certain balance among the 
two, possibly depending on the feedback mechanisms perform¬ 
ing differently at different luminosity and mass regimes. 

Deeper data with HST/WFC3 UV channels would be needed 
to examine the spatial distribution of the ionizing radiation (as a 
probe of a stellar ionizing emission), and/or deeper spectroscopic 
observations would be necessary to study the covering fraction 
of neutral gas through absorption fine analysis (e.g., Heck man et 
|al. 12011[|Borthakur et al. [2014| . 


4.3. LyC candidate: Ion2 

Another_source at z = 3.218 (CANDELS ID 18320, [Guo et| 
|al. ||2013| with a significant excess in the ultraviolet spectrum 
has been identified in the same field with a U-band magnitude 
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U = 26.85 ±0.06 and color (J/-fi)=1.66±O.O70At this redshift 
the U -band probes the wavelength range A < 948A rest-frame, 
and any significant signal suggests a LyC leakage. From the ob¬ 
served (U - #) color and following the methodology described 
above, we derive P(fesc=0)=0.0% (Fig. [5]), which confirms that 
this source is a LyC candidate (named here Ion2). This result 
also supports the non-zero emitting flux blueward of the Lyman 
limit that was tentatively detected in the VIMOS low-resolution 
spectrum, R = 180 (see Fig. but with a low S /N <2 ( |Balestra| 
jet al. 120101 ). 


By examining the morphology of the system in more de¬ 
tail, two distinct objects are visible at a separation of 0.2" in 
the HST/ACS imaging (corresponding to a physical separation 
of 1.5kpc if they are at the same redshift, Fig. [5]), both of them 
showing a drop blueward of the V 606 -band. 

In the following, we refer with a) to the brighter and with b) 
to the fainter component. They have been detected as a single 
object both in the optical ACS catal og v2.0 (Giaval isco et al.| 
H |2004| and see |Vanzella et al. ~||2009 ) and in the near-infrared 
CANDELS catalog ( Guo et al. ||2013| ). They are not detected 
in the X-ray 4Ms Chandra observations, corresponding to lcr 
limit of Lx < 4.1 x 10 42 erg s~ l . That these are two sources 
further complicates the interpretation. In the following sections 
we investigate the nature of the system in more detail. 


4.4. Photometric and spectroscopic confirmation of LyC 
candidates 


4.4.1. LyC emission 


Another useful feature is detected in the medium-resolution 
spectrum obtained with VIMOS (R = 580) that was published 
in |Balestra et al. | ( [2010| ). The Lya line shows a clear and atypi¬ 
cal asymmetric blue tail, at variance with the typically opposite 
shape toward red wavelengths (|Shapley et al. |2003| |Verhamme] 


et al. |2008| j Vanzella et al. 2009[ ). We performed a dedicated 

reduction of the VIMOS quadrant containing Ion2 following the 
reduction prescription described in [Vanze l la et al. [{2014b| ). The 
two-dimensional spectrum is shown in Fig.[ 6 ](panel d), in which 
the Ciii]/U909 emission line at z - 3.212 is clearly identified 
(S/N= 8 ), and the Lya was detected at the first and second order 
(the latter with a double spectral resolution, R ^ 1200). In par¬ 
ticular, the second-order detection shows a clear double-peaked 
emission with dv - 550 ± 90 km s~ x (see Fig. [bj inset of panel 
b). There are at least two possible explanations for the two lines 
in emission around the Lya position: (1) the two sources iden¬ 
tified in the ACS images are at a very similar redshift and both 
with Lya in emission and detected in the MR spectrum; (2) one 
of them has an intrinsic double-peaked Lya line in emission (as 
observed, e.g., in |Karman et al. |20l4} |Kulas et al. |20~l2] ). It is 
worth noting that Ly/3 and Lyy lines have been clearly detected 
in absorption in the LRb spectrum; these detections are possible 
only for the brighter of the two components. The fainter compo¬ 
nent has an estimated magnitude at the wavelength interval en¬ 
compassing the Ly/3 and Lyy absorption lines (i.e., - # 435 -band) 
of L 435 = 27.25 ± 0.24, which is too faint to allow a detection 
of the continuum in the VIMOS LR spectrum with only four 
hours of integration time (expected to have S /N ~ 1, running 
the ES O/Exposure Time Calculator^] and therefore the detec¬ 
tion of any absorption line is prohibitive, too. This implies that 
they are associated with the brighter component (at ZLya = 3.21), 


4 The #-band data used here come from the ultradeep VIMOS obser¬ 
vation of GOODS-S, Nonino et al. in preparation. 

5 http://www. eso. org/observing/etc/ 


that is, two magnitudes brighter in the # 435 -band. The uncer¬ 
tainty remains whether the (fainter) source (b) is at the same 
redshift as the main object or not. If it is at a very similar red¬ 
shift, then it would be interesting to also investigate its ionizing 
emission. Conversely, a lower redshift source z « 3.2 would 
be less interesting and could contribute to the observed U -band 
magnitude (but not to the LR spectrum, being too faint). The 
HST/ACS resolution allows us to perform dedicated photometry 
of both components, in particular, components a) and b) were 
subtracted from each ACS band (see Sect. 4.6), and the pho¬ 
tometry was properly derived by running SExtractor. The two 
components show very similar photometric behaviors; we derive 
the following ACS BViz magnitudes for the (brighter) source a): 

25.34 ± 0.13, 24.72 ±0.05, 24.66 ± 0.02, and 24.64 ±0.12, and for 
the fainter source b): 27.25+0.24, 26.63+0.09, 26.41+0.04, and 

26.35 + 0.20. Source b) also shows a clear drop in the # 435 -band, 
appreciable in Fig. [5] If interpreted as due to the IGM decre¬ 
ment, it would be compatible with the redshift and the decre¬ 
ment also observed in the a) companion (ZLya = 3.218). The 
(#435 - V606) colors of the two components are practically the 
same, 0.62+0.14 for a) and 0.62+0.26 for b), while the # 435-/775 
drop is slightly larger for the fainter one, supporting its high- 
z nature. The other possible interpretation for the drop is the 
Baimer break at z - 0.2. However, at this redshift and given 
the spectral coverage in the range 3500 - 10000A, typical fea¬ 
tures such as the [Oii]T3727, [Om]/I4959 - 5007, H(3 , and Ha 
would be easily detectable in the spectra if it were an emission 
line object. Regardless of the redshift of the fainter component, 
the estimated 17-band magnitude of 26.85 + 0.06 for the system 
cannot be explained by the b) source alone, even assuming that 
the b) component maintains the same #435 = 27.25 + 0.24 mag¬ 
nitude in the 77-band as well. 

The two distinct and close objects, the two lines in emission, 
and the drop blueward of the V 606 -band in both objects support 
the interpretation that the two sources are at very similar redshift 
with a velocity difference of Av ^ 550 + 90 km s~ l , with the 
two Lya lines resolved in the R = 1200 spectrum (they are not 
spatially resolved with the given slit orientation). Under this as¬ 
sumption, the system remains a convincing LyC candidate, with 
P(fesc=0)= 0. 


4.4.2. Low gas-covering fraction ? 

|Heckman et al. | ( |2011| ) suggested a link between the low- 
ionization [Cn]/U334 (and [Sin]/il260) absorption line and the 
transparency of the interstellar medium to ionizing photons, such 
that a non-zero residual flux in the [C n]/U334 line is a necessary 
condition to have a low optical depth at the Lyman edge (see 
also |Jaskot & Qey |2014| ). |Borthakur et al. | ( |2014| ) found similar 
features in a local Lyman-break analog with Lyman-continuum 
leakage. In the present case, the absence of the low-ionization 
absorption lines of silicon and carbon (see Fig. [ 6 ]), especially 
the faintness of the multiple atomic transitions (i.e., the evident 
non-zero residual flux) of silicon [Sin]/11260,1304,1526, and 
[C n] /U3 34 suggests a low covering fraction of neutral gas (see 
Jones et al. 2013} . The available spectral resolution prevents us 
from performing a more quantitative study, but the absence of 
these lines supports a possible high transmission along the line 
of sight, consistent with the potentially high escaping ionizing 
radiation we inferred indirectly from the MC analysis. 

As concluded in the previous section, it is plausible that 
the two components are at the same redshift, both with Lya in 
emission. Another piece of information comes from the semi- 
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Fig. 6. Low- (R = 180) and medium- (.R = 580) resolution VIMOS 
spectra of Ion2 in panels a) and b). The two spectra overlap in the 
range (4800-6800A). Several typical low- and high-ionization atomic 
transitions are reported. The Lya emission in both spectra and the Ly/3 
and Lyy absorption features in the low-resolution spectrum (panel a) 
are clear. The Lyman continuum blueward of the Lyman limit (LL) is 
indicated in panel a). Panel b) (bottom left) shows the zoomed one¬ 
dimensional Lya line extracted from the medium-resolution spectrum 
(R=580). Its atypical blue asymmetric tail is clear. The insets in panel b) 
show the two-dimensional Lya line at first (top) and second (bottom) or¬ 
der. The latter, with double spectral resolution, R ^ 1200, clearly shows 
the double-peaked structure of the line. In the top right panel (panel c) a 
zoom-in of the position of C m]/tl909 is shown at z = 3.212. The insets 
in the panel d) show the region in which the Cn]/11334 absorption line 
would lie: absorptions consistent with noise fluctuation are present (the 
gray stripe indicates the position of the sky emission line, 5511 A). 


forbidden doublet (here unresolved) nebular C iii]/LU906- 1909. 
In particular, the Zcni] = 3.212 and the redshift of the bluer 
peak of the Lya structure, ZLya,biue = 3.211, are fully compat¬ 
ible within an uncertainty of 10km s~ l , while the red peak is 
redshifted by about 550 km s~ l (Ziyajed - 3.218 ± 0.001). If we 
adopt the Cm] transition as a proxy of the systemic redshift, 
then the emerging Lya photons at the systemic redshift could 
further support a scenario of a low neutral gas column density 
that would be consistent with the possible ionizing leakage (the 
resonant scattering is reduced and the Lya photons escape not 
far from the resonance frequency, close to the systemic velocity, 
see |Verhamme et al. |2014t|Schaerer et al. |2011| ). If this is the 
case, the escaping ionizing and Lya photons could share a sim¬ 
ilar (or have the same) physical path through the cavities in the 
interstellar medium ( [Behrens et al. |2014| . 

4.5. Signature of variability: faint AGN component ? 

The ACS /5435 and VLT/VIMOS 5-band observations have been 
obtained in two epochs, in 2002 and 2010, respectively. The 
two filters encompass a very similar wavelength intervals with 
a small difference in their shape (see Figs. [4] and [5]). We quan¬ 
tified the magnitude difference (B V imos~B 435) to be in the range 
0-0.1 (and average +0.03) by convolving the two filters with a 
large control sample made of BC galaxy templates. The differ¬ 
ence of the observed magnitudes A B=(B V imos ~ #435) at the two 
epochs (corrected for the average offset) and the uncertainty cr^ B 
can be computed. The quantity BVAR=ABS(AB / ct^b) provides 
an estimate of the significance of the time variation between the 
two epochs (e.g., [Villforth et al. |2010| . The rest-frame interval 
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Fig. 7. Example of the GALFIT fitting in the / 775 -band of the two com¬ 
ponents of Ion2 . In the top and bottom panels the a) and b) compo¬ 
nents have been subtracted. The best solution provides an effective ra¬ 
dius R e ^ 1.1 pixels. 


of time probed at redshift ~3 is ^ two years. As a sanity check, 
we calculated this quantity for the known AGNs in the GOODS- 
S field ( [Szokoly et al. |2004]|Xue et al. |2011| ), and for several of 
them the BVAR ranges between 6-10, meaning that the variation 
is evident and significant at more than six sigma (as expected for 
AGNs, especially in the ultraviolet, e.g.,[Cristiani et al.] T 997j ). 
This will be explored in a forthcoming work on the full photo¬ 
metric sample in the GOODS-S field. Here we limit our analy¬ 
sis to the two examples reported above. The derived BVAR for 
lonl and Ion2 is 1.3 and 3.0, respectively (corresponding to a 
magnitude difference A B = +0.26 and A B = +0.42). The two 
sources are isolated and not subject to photometric contamina¬ 
tion from other nearby objects. While for lonl there is no signif¬ 
icant evidence for variability, in the case of Ion2 it is significant 
at a 3 <x level. Since the variability is a sufficient condition to 
confirm nuclear activity, its detection in the case of lon2 sug¬ 
gests that an AGN component is highly likely. We recall that 
the two sources are currently not identified with deep 4Ms X- 
ray observations (X ue et a l. [2011[ ), corresponding to l<x limit 
of L x < 4 x 10 42 erg s~ L for both. In addition, the typical high- 
ionization spectral features of nuclear activity (e.g., N V T1240, 
Civ T1550) are absent as well. However, in a fainter luminosity 
domain the situation could be different. [Silverman et al. 1 |2010| ) 
have found X-ray selected AGNs at 2 < z < 3.6 with optical 
spectra fully compatible with those of typical star-forming galax¬ 
ies, that is, without high ionization emission fines and/or without 
any fine in emission (see also Civano et al. 2011) . Without the 
X-ray imaging the AGN component in these systems would not 
have been identified. The UV slopes of these sources is also com¬ 
patible with the UV slope observed here. [Vanzella et al.[(j2010b| 
also reported on two AGNs at z — 3.462 and z = 3.466 ( palestra 


|et al. [ 2010 | ) with ionizing radiation directly detected in the U- 
band ultradeep VIMOS observations, both showing a blue color 
( 1*775 - z Z 85o) ^ 0. The two sources lonl and lon2 are ~ 1 magni¬ 
tude fainter, not detectable in the deep X-ray observations, and 
may represent a regime in which the stellar and non-stellar ultra¬ 
violet emissions are similar. 


4.6. Spatially resolved emission 

High spatial resolution ACS imaging has been fitted with GAL¬ 
FIT ( [Peng et al. |2010| to extract the basic PSF-corrected mor- 
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phological parameters. The assumed PSF in the HST images has 
been derived empirically from observed stars close to the targets. 

First, we note that lonl and Ion2 (both the a) and b) com¬ 
ponents for the latter) show a compact uncomplicated shape, 
so a good fit is reached by adopting a simple Gaussian profile 
(Sersic model with n = 0.5) and leaving the effective radius Re , 
the axis ratio B/A, the coordinates X,Y, the magnitude, and the 
position angle (PA) as free parameters (very similar results are 
also obtained by leaving the Sersic index n as a free parameter). 
B/A and PA are not relevant parameters in the fit, since all com¬ 
ponents have a circular symmetry, therefore we fixed them to 
B/A=l and PA=0. The coordinates of the sources and their mag¬ 
nitudes were fixed to the best guesses derived from SExtractor. 
Then GALFIT was run many times by varying the Re parame¬ 
ter from 0.025 to 2.0 pixels. An example is shown in Fig. [7] in 
which the two components of Ion2 were modeled and subtracted 


from the / 775 -band (^ 1800A rest-frame). The behavior of the 
residuals at the position of the source is shown as a function of 
the effective radius (right panel of Fig. [7]). The minimum of the 
residuals for Re - 1.1 - 1.2 pixels is clear (lpix=0.03"), and the 
residuals map provided by GALFIT, as a result subtracting the 
best-fit model from the original galaxy, do not show significant 
structures (middle panels of the same figure). 

We performed this exercise in all the ACS bands (^ 1000 - 
2000A rest-frame), both for lonl and Ion2. Both sources (and 
the two subcomponents of Ion2) are spatially resolved in the ul¬ 
traviolet, with effective radii of 200 and 300 parsecs (± 70 pc), 
respectively. While in the previous section we discussed the pos¬ 
sible presence of an AGN component, especially from the vari¬ 
ability of Ion2 , the spatially resolved emission suggests that stel¬ 
lar UV light also contributes to the emerging spectrum. Unfortu¬ 
nately, high spatial resolution imaging in the ionizing domain is 
not available at the moment. Moreover, the compactness of the 
central star-forming region within a few hundred parsecs and the 
leakage of ionizing photons has recently been reported in a local 
starburst galaxy ( [Borthakur et al. |2014| ), in which a high volume 
density of young stars can generate intense winds able to propel 
the surrounding interstellar neutral and ionized gas outward at 
velocities in excess of 1000 kms~ l (Heckman et al. 20lTJ». The 
two sources discussed here could be the higher redshift counter¬ 
parts and will merit future investigation of the possible presence 
of winds. 


5. Discussion and conclusions 

The rest-frame far UV colors of star-forming galaxies, which at 
z > 2.7 are observed in the optical window, are powerful diag¬ 
nostics of the forming stellar populations and of the ISM. They 
are also very useful tools for selecting the galaxies themselves 
by means of dedicated color-color diagrams, such as the Lyman- 
break technique ( Giavalisc o |2002| ), or by means of photometric 
redshift selection. The spectral region most often targeted for 
the construction of colors typically straddles the 912A Lyman- 
continuum discontinuity and the Lya resonant line at 1216A. In 
general, the dispersion of the observed colors results from the 
contribution of a number of independent factors, which include 
the following: 

1) the intrinsic scatter of the UV SED of galaxies. This con- 
trbution primarily depends on the scatter of the integrated colors 
of massive stars and on the scatter of reddening by dust. The 
former term is dominated by variations of the age of the star- 
burst and variations of the IMF, the latter by variations of the to¬ 
tal amount of dust (which is commonly parameterized in terms 


of the E(B-V) color excess) and of variations of the reddening 
curve; 

2 ) the scatter of the cosmic opacity, which consists of the 
Lyman-alpha forest and of DLAs and LLSs. Detailed analyses 
have been performed in this work by running MC simulations 
based on the prescriptions of 108 and 114; 

3) the scatter of the amount of escaping ionizing radiation, 
which depends on the HI column density of the ISM of the 
source along the line of sight. There are no firm constraints on 
this term, but its effects on the observed color are certainly weak 
because the number of leaking ionizing photons is small in abso¬ 
lute terms compared to non-ionizing photons. Another cause is 
the scatter of broad-band colors constructed with passbands that 
sample the spectral regions to the blue and the red of the Lyman 
limit are much more strongly affected by small differences in the 
redshift of the targeted galaxies, which determine the amount of 
non-ionizing flux that enters the blue passband. This amount of 
flux is generally much larger than the potential leaking ioniz¬ 
ing flux, and so is the corresponding contribution to the scatter 
of colors. Nevertheless, knowledge of the exact position of the 
Lyman limit (spectroscopic redshift) has allowed us to define a 
method for identifying the signature of a possible escaping ion¬ 
izing radiation; 

4) photometric scatter. This term depends on the depth of the 
photometry used to construct the colors, and thus on the experi¬ 
mental configuration. It is usually accurately measured as a func¬ 
tion of the luminosity of the source by means of Monte Carlo 
simulations. It can provide a comparatively strong contribution 
to the scatter, especially in colors that include spectral regions 
where the source is faint, for instance, blueward of the Lyman 
limit, the Lyman-alpha forest, or both; 

5) confusion scatter. This term arises from photometric con¬ 
tamination by nearby spuri ous interlopers at low er redshift than 
the targeted source (e.g.,|Vanzella et al. 201 ()a| ). It depends on 
the angular resolution of the images used to construct the col¬ 
ors, and thus on the instrumental configuration and ground-based 
seeing (if applicable). Even if the interloper is a very faint galaxy, 
the contamination can be severe in colors that include spectral 
regions where the source is faint, for example, blueward of the 
Lyman limit, the Lyman-alpha forest, or both. 

We analyzed the effect of an escaping ionizing radiation on 
color selection criteria by running dedicated Monte Carlo simu¬ 
lations. On the basis of similar simulations anchored to the ob¬ 
servations, a method of selecting LyC continuum emitter candi¬ 
dates from multiband photometry was also discussed and applied 
to a sample of galaxies. 

These are the main results: 


1. The redshift evolution of the mean free path of ionizing pho¬ 
tons through the intergalactic medium and a proper treat¬ 
ment of the ionizing continuum (galaxy template) indicate 
that there is no significant effect of LyC leakage on classical 
broad-band photometric selection techniques (e.g.,U n GR or 
BVI ), at variance with what has been reported by|Cooke et 
|aT p0l4) . The effect of an ionizing emission is more sig¬ 
nificant as the width of the band encompassing the Lyman 
edge decreases and approaches the expected mean free path 
at the given redshift (being maximized in the narrow-band 
surveys or deep spe ctroscopy, e.g.,|Mostardi et al. |2013||Gi- 
|allongo et al. |2002{|Shapley et al. |2006| ). Under the assump- 

tion that the IGM prescription of 114 is a good representation 
of the IGM structure, no significant deviation was found in 
the non-ionizing part of the IGM attenuation (912 - 1216A) 
either, that is, the scatter is not larger than 0.4 mags. Over- 
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all, the median color tracks do not differ significantly if fesc 
= 0 or 1, even in the most extreme stellar ionizing emitters 
(age= lMyr, E(B-V)=0, Z/Z Q = 0.02). 

2. A method for photometrically selecting LyC-emitter candi¬ 
dates was described and applied to a sample of 35 galaxies 
belonging to the GOODS-S field. Two candidate LyC emit¬ 
ters were identified: one known LyC emitter at z = 3.795 was 
successfully recovered ( lonl ). Another source with possible 
LyC leakage was reported (Jon2). The probability of a null 
escaping ionizing radiation is low (< 5%) for both sources. 
Detailed analysis based on morphology, variability, and the 
ultraviolet excess suggests that a probable contribution from 
faint AGN activity could dominate the LyC emission and/or 
have played a crucial role in making the medium transparent. 
The spatially resolved stellar UV emission also suggest there 
could be a contribution by stellar radiation. 


There are caveats for the method described in this work. 
Lirst, the probability P(fesc=0) depends on the fraction of the 
band used to probe the LyC region. A large fesc could be present, 
for example, but might be missed because only a small fraction 
of the band probes it. This effect can be attenuated by includ¬ 
ing bluer overlapping bands that ensure a continuous redshift 
coverage (as is the case in the typical multifrequency surveys 
such as GOODS/CANDELS). Second, the method tends to se¬ 
lect relatively strong LyC emitters on average (it depends on the 
width of the filters used and the presence of LLSs/DLAs IGM 
along the line of sight) because they are the most effective in the 
broad-band contribution. Third, the photometric accuracy is cru¬ 
cial in these analyses to identify the position in the color-color 
plane with good detail and determine P(fesc). The higher the 
photometric S/N ratio in the relevant bands, the more accurate 
the identification of LyC candidates with this method. The inter¬ 
play among these limitations is included in the performed MC 
simulations, however. 

Despite these limitations, we have demonstrated that regard¬ 
less of the nature of the source (AGNs, LyC-galaxy, spurious 
detection, etc.), this type of objects can in principle be automat¬ 
ically identified with this technique. The color excess in the UV 
coupled with morphological information has previously been 
used to select QSOs/AGNs (e.g. , |Lontanot et al. |[2007| |Casey| 
|et al. |2008| |Salvato et al. |2Q09) and is conceptually similar to 
what we discussed here. However, the analyses performed here 
are based on the detailed IGM stochasticity and focused mainly 
on the LyC part of the spectrum. 

Lrom point (1) above it follows that the broad-band photo¬ 
metric selection of high-redshift galaxies does not prevent se¬ 
lecting galaxies with LyC leakage. Nonetheless, a clear sample 
of ionizing galaxies has not been identified yet. One reason could 
be that the stellar ionizers exist in a much fainter and still unex¬ 
plored luminosity domain, while relatively bright (L > 0.5L*) 
ionizing sources are rare and require efficient physical mech¬ 
anisms to increase the transparency of the interstellar medium 
and favor ionizing photons to escape. The variability detected in 
Ion2 and the difficulty in reproducing (with galaxy templates) the 


observed U-band magnitude of lonl (A < 830A) given its ob¬ 
served ultraviolet slope suggest that the two sources could host 
a faint AGN that possibly contributes to the LyC emission. X- 
ray detected ( L x > 10 42,5 erg s~ x ) AGNs with blue UV-slope 
and galaxy-like spectra have been reported in the literature (e.g., 
|Civano et al. |201l||Silverman et al. |2010| ), and the two s ources 
reported here could be AGNs at fainter luminosities (e.g., Maoz 
2007| ). In other words, while the more massive systems would 
need an efficient feedback (possibly provided by an AGN) to 
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clean the medium and decrease the column density of the neu¬ 
tral gas, the fainter and still unexplored sources could regulate 
the opacity of the interstellar medium with supernova feedback 
only if they indeed contain a large portion of the LyC emitter 
( [Wise et al. |20l4l|Kimm & Cen |2014) . 

In addition to the possible presence of a faint AGN com¬ 
ponent, the spatially resolved emission in the ultraviolet 1000- 
2000A rest-frame suggests that the stellar emission contributes 
in this wavelength domain. It is beyond the aim of this work 
to decompose the two components, we only note that the two 
might also coexist in the LyC domain. In particular, stellar ion¬ 
izing photons could share a very similar path of the interstellar 
and circum galactic media of these photons arising from the nu¬ 
cleus, making the system a sort of hybrid-ionizer (stellar and 
non-stellar). In this respect, high spatial resolution imaging with 
HST (WLC3 UV channels) of confirmed LyC emitters would 
provide precious information about the spatial distribution and 
the relative balance of the two sources of ionizing radiation. 

If there were a faint AGN component in this type of sources, 
it might also explain the recently observed extreme ultravio¬ 
let flux ratios between ionizing and non-ionizing wavelengths 
in various Lyman alpha emitters reported in the literature (e.g., 
/ 1500// 900 < 1-2), e.g., |Nestor et al. | ( |2013| ), [Mostardi et al 


( |2013| ). Lor example, Ion2 is a Lyman alpha emitter without any 
typical feature, which suggests that it is an AGN (from X-ray 
and spectral features), but its variability suggests that there is an 
AGN component. We therefore remark that it is important to in¬ 
vestigate the AGN nature at faint luminosity limits in this type 
of studies. 

Linally, convolving the stochasticity of the IGM with galaxy 
templates anchored to the observed UV slopes and comparing 
this with the observed colors proved to be an efficient automatic 
tool to determine the reliability of spectroscopic redshift mea¬ 
surements or to identify photometric contamination from lower- 
z sources in the UV. In these cases, an addiitonal flux in the ultra¬ 
violet is present and is therefore recovered as a possible ionizing 
leakage. This is particularly useful for large spectroscopic red¬ 
shift surveys coupled with deep and high spatial resolution imag¬ 
ing, such as VUDS (Le Levre et a l. |2Q15| ) and/or the VANDELS 
survey (ESO public survey^]), which cover the CANDELS fields. 
We defer the application to these large datasets to a forthcoming 
work. 
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